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Abstract 

In this paper we present two numerical schemes of approximating so- 
lutions of backward doubly stochastic differential equations (BDSDEs for 
short). We give a method to discretize a BDSDE. And we also give the 
proof of the convergence of these two kinds of solutions for BDSDEs re- 
spectively. We give a sample of computation of BDSDEs. 
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1 Introduction 

Since Pardoux and Peng introduced backward stochastic differential equation 
(BSDE), the theory of which has been widely used and developed, mainly be- 
cause of a large part of problems in mathematical finance can be treated as a 
BSDE. However it is known that only a limited number of BSDE can be solved 
explicitly. To develop numerical method and numerical algorithm is very help- 
ful, theoretically and practically. Recently many different types of discretization 
of BSDE and the related numerical analysis were introduced. 

On the other hand, Paroux and Peng Q introduced a new class of backward 
stochastic differential equations-backwarddoubly stochastic differential equa- 
tions and also showed the existence and uniqueness of the solution of BDSDE. 
But until now little work is devoted to the numerical method and the related 
numerical analysis. Here following the approach of Mcmin, Peng and Xu 
we present two numerical schemes of approximating solutions of BDSDE, and 
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proved the convergence of these two kinds of solutions for BDSDEs, respectively. 
First of the proofs makes use of and extends Donsker-Type theorem. 

This paper is organized as follows. In section 2, we introduce some fun- 
damental knowledge and assumptions of BDSDEs. In section 3, the discrete 
BDSDE and solutions are presented. In section 4. wc will give our main results: 
the proof of convergence of numerical solutions for BDSDEs in two different 
schemes. 

2 Some Preliminaries 

Let (r2, P) be a complete probability space, and T > be fixed throughout 
this paper. Let {Wt, < t < T} a,nd {Bt, < t < T} be two mutually 
independent standard Brownian motion processes, with values respectively in 
M-^and in M', define on (f7,J^,P). For each t e [0,T], we define 

Tt = T^' V Tt^T 

where for any process {rjt}, J-l'^ ~ <j{rir -~ rjs] s < r < <}, ~ J'q t- 

For any n G N, let M'^{0,T;W"-) denote the set of (classes of dP x dt a.e. 
equal) n dimensional jointly measurable random processes {ipt; t G [0, T]} which 
satisfy: 

T 

(i) . E / \ipt\^dt < oo 



(ii) . (fit is J^t-measurable, for a.e.t G [0,T]. 

We denote similarly by iS^([0, T]; R") the set of continuous n dimensional 
random processes which satisfy: 

(i) . E(supo<4<r \ftf) < oo 

(ii) . ipt is J-j-measurable, for any t e [0,T]. 
Let 

/ : f7 X [0, T] X R'^ X R'''"* ^ R'= 
g ■.nx[0,T]xR'' X R*^^'^ ^ R'=^' 
be jointly measurable and such that for any {y, z) £ M.'^ x R'^^'', 

f{;y,z)eM\0,T;R'') 

g{-,y,z)eM\0,T;M.''>^') 

We assume moreover that there exist constants K > and < a < 1 such 
that for any {uj,t) € 17 x [0, T], (j/i, zi), (j/2, ^2) G R'' x R'^^'', 

(HA) \f{t, y,,z,) - fit, y2,Z2)\ < Ki\y, - y^l + \\z, - Z2II) 
\\g{t,yi,zi) - g(t,y2,Z2)\\ < K \yi - yal +a\\zi - Z2II 
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Given ^ G L^(r2, P; K^), we consider the following backward doubly 
stochastic differential equation: 

T T T 

Yt=i + j f{s, Ys, Zs)ds + j g{s, Ys, Zs)dB, - j ZsdWs, 0<t<T. 
t t t 

We note that the integral with respect to {Bt} is a "backward Ito integral" 
and the integral with respect to {Wt} is a standard forward forward Ito integral, 
see Nualart and Pardoux Q • 

Here we mainly study the case when Brownian motion is one-dimensional. 
Now we consider the following 1-dimcnsional BDSDE 

T T T 

Yt^^ + J f{s,Ys,Zs)ds + J g{s,Ys,Z,)dB,- J Z^dWs (1) 
t t t 

and the terminal condition is i/t ~ £, = $(Wt), where $(•) is a functional of 
Brownian motion {{Bg, Ws)o<s<t} , such that G L'^{Tt)- Particularly, if /(•), 
g{-) are not relative to t, ([1]) changes into: 



T 



Yt=^ + j f{Y,, Zs)ds + J g{Ys, Z,)dB, - j Z,dW, (2) 



t 



3 Numerical Scheme of Standard BDSDE 
3.1 The Structure of Numerical Solution 

When n e N is big enough, we divide the time interval [0,r] into n parts: 
= io < < • • • < tri = r, (5 := tj - tj_i = ^, for l<j<n. 

Now we define the scaled random walk {-B"}, {VK."} , by setting Bq = = 

0, 

[t/&] [t/S] 

Bt = wr = vsY,P7, o<t<T 

where {£"}"^i, are two mutually independent Bernoulh sequences, which 

are i.i.d. random variable satisfying 



= 0.5 
0.5 



Obviously, i?", W" are both JF^ -measurable processes who take discrete val- 
ues, denote B^ = Bl\,Wf = Wl], we get = £m, W^" = E 

m—l r—1 

And we define the discrete filtrations fj" = a{(3i, (3j} = (t{Wj",0 < t < 
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= <T{/3i,...,/?,}Va{£,+i,. ..,£„} = a{W,",0 < t < i,} V ^{5,", < 

t<T}, 

Then, on the small interval [tj,tj+i], the equation 

tj+i tj+i tj+i 

Yt^=Yt^^,+ J f{s,Ys,Zs)ds+ j g{s,Ys,Zs)dB,~ j Z,dW, (3) 

can be approximated by the discrete equation 
i.e. 

y7 ^ y?+i + fi^^^y?^ + 5(^^+1, y^+i, ^;-Vi)V^£,+i - (5) 

For sake of simplicity, here we just consider the situation in which /, g are 
not relative to t. 

Lemma 3.1. Let be a given f/jYij+i "measurable random variable. Then, 
when 5 < 1/k, there exists a unique C/^-measurable pair {y'^,z^) satisfying the 
equation: 

yl = + /(y"' ^1)^ + 5(2/"+i, - (6) 

Proof. We set Y+,^ = y^^, |ft+,=i,F^+i = y^+i |/3,+,=-i, 2/+ = y," |e,+,=i 

) Uj — Vj |e3 + i = -l • 

Both Y+^„Yr^, are ^^j^-mcasurablc. Then equation ([6]) is equivalent 
to the following algebraic equations: 

yt = Y+, + f{yt.-4)5 + 3(1^+1, - 

yt = ^+1 + /(y/' + 9{Y-,„z-^,)V~5 + z+VS 
yj = Y+^^ + f{yj.zj)5 - g{Y+^„z+^,)VS- zjVS 

yj = Yf+^ + f{yj.zj)5 - 1, z-^,)V5 + zjVs 
Solving these equations, we can get 
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That is to say: 

4 = - ^.+1) + ^[.9(^+1' 4+1) - g{Y-,i,z7+i)]e,+i 

y? - fiy?' ^7)^ = + V+i) + ^l9iY;+„z++^) + 9{Y-,,z-_^,)]e,+, 

We can simulate a sample path of {sj}, then we calculate the corresponding 
BSDE along with the sequence. It is indeed a kind of Monte-Carlo method. 

Example 3.1. If f{y, z) ~ ay + bz, 

J l-a5 



The calculation begins at the terminal time f„ = T, with yJJ = which 
is given and the problem is how to determine Z„. Here we choose the way of 
setting Zt = Vlr, i-C. Zt = dxYx- 

Example 3.2. For simphcity, we suppose a linear type, f{y,z) = 0, g{y,z) = 
ay + bz, Yt = Wt, then Zt = VFt- 
We have 

= Y+ + (aY+ + bZ+)^r5 - Z^^^VS 
Y+_, = Y+ + {aY+ + bZ+)V5 + Z+_,V5 
Y--1 = Y+ + {aY+ + bZ+)^r5 - Z-_^V5 
Y-^r = 1;+ + (ar+ + + Z-_^^r& 



I.e. 
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Z--, = ^""^^^"" (1 - aV&) + \{Z+ Z-) 

After Zn is calculated, Yj and can be backwardly step by step, following 
the way mentioned above. 

On the other hand, taking conditional expectation on ([51, it follows that 

At the terminal time t„ = T , consider the mapping '^{z) — z — ^[g{Yj^-, z) — 
g(y_, z)]£j-l-i, from the property of g, we obtain that the derivative of ^(z) 
on z is 1, which implies that the mapping 'I'(z) is a monotonic mapping. So 
there exists a unique value Zn-i s.t. Zn-i = ^~^(^^ ^ ^ ) + ^[.^(^^i z„_i) — 
Zn-i)] holds. Consider the mapping Q{y) = y — f{y, Zj-)S from the Lips- 
chitz property of /, we obtain 

(e(y)-e(2/'),2/-y') > {I-SK) \y-y'f>0 
which implies that the mapping Q{y) is a monotonic mapping. So there exists 
a unique value y s.t. e(y) = E[y;Vi I S^j] + VSE[g{y^_,,, z^+^) \ g])]ej+i holds, 
i.e. = Q-^nv'j+i I S",] + ^/^E[g(y7+i,z7+i) | Q],]e,+i). 

Remark a. The existence of the solution of discrete BDSDE only depends 
on the Lipschitz condition of / on y. In fact, if / does not depend on y, we can 
easily get Q~^{y) ^ y + f{zJ)S, And very obviously, if g does not depend on z, 
{z"} can be also easily got. 

Remark b. In general, if / nonlinearly depends on y, then 8(y) can not be 
solved explicitly, so sometimes we can introduce the following scheme, we set 
Yrp =27" = and starting from j = n ~ I, backwardly solve 

y7 = y5Vi + /(E[y5Vi l^:,]^^)'5 + 5(y;Vi,^i)Vfc,+i-z5'V^/3,+i (7) 

or equivalcntly, 

z] = ^E[y]+J3,+, I g^^] + E[5(y;Vi: I 5J,]e,+, (8) 

y] = E[y^i I g^^] + fmy]+i I Sf,l^l)S + E[.g(y5Vi, (9) 
to approximate the solution of 6(y) = E[?/"^^ | C/j^] + \/SE[g{yJ_^i, Zj-_^_i) \ 

g-]e,+^. 
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3.2 Monte-Carlo Method 



For Forward-Backward SDEs, 



dXs = b{s,X,)ds + a{s,Xs)dWs. 
Xs = x, 0<s<t. 



(10) 



dn = f{s, Z,,) ds - dM^„ 



(11) 



where 6, a, f and $ satisfy usual assumption, by Theorem 4.1 of [3 There exist 
two function u{t,x) and d{t,x), such that the solution (F*'^, Z*'^) of BSDE is 

r,*'^ =u(s,X*'^), Z*'^ = cr(s,X*'^)d(s,X*'^), < < s < T, dP (8) ds a.s. 

The solution of the BSDE is said to be Markovian. So it is naturally to solve 
the equation based on a binomial tree of Xg. 

Example 3.2. If Xs = Ws, the solution of BSDE is 

= u{s, Ws), Zs = d{s, Ws), i < s < T, dP (g) ds a.s. 

As for BDSDE, the structure of BDSDE is different from BSDE, that the 
solution is not generally in the form of Yt = Wt, Bt—Bi), Zt = ip{t, Wt, Bt~ 
Bf ), even though / and g are deterministic functions. 

Example 3.3. 



Therefore, Yt is path dependent on B., So it's impossible to solve the solution 
on the nodes of the coupled binomial trees. 

If we simulate a sample path of Bt , it becomes a classical numerical scheme of 
BSDE follow the path, which is indeed a Monte-Carlo method, and the solution 
surface will vibrate with the sample path of Bt . Yang [221 gives some comparison 
examples of numerical solutions and explicit solutions. □ 

^The reason we don't include the examples in this paper is that arXiv rejeet figures of large 
size. 



dYt 

Yt 



tdBt-ZtdWt 



0. 



t, 



The solution is 
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3.3 Associated SPDE 

For each {t, x) e R+ x R^, let {X*'^, i < s < T} be the sohition of the SDE: 

S S 

X*'^ = x + J h{Xl-'')dr + j <7{X*-''=)dWr, t<s<T. 



The foUowing BDSDE 

T T T 

F^*.- = h J fiX^-'^X'"'^ Zl.-'=)dr+ y .9(X;^^ y;■^ Zl'^)dBr-J t<s<T. 

s s s 

Under assumption (HI) has a unique solution {Y^-^ , Z^-'^), and under some 
suitable conditions, 

u{t, x) = y/'"', (i, x) e [0, T] X R'' 
is the unique solution of the following SPDE: < i < T, 

u{t, x) = u(T, x)+ I [Lu(s, x)+f{x, u{s, x), (crVw)(s, x))] ds— / 17(2;, u{s, x), (ctVu)(s, a;)) dBg 



Note that u{t,x) depends on B{-) indeed. 
Example 3.4. / = 0, g = 1, 



u{t,x) = u{T,x) + / Lu{s,x)ds+ / dBs- 



It Jt 
Wt itself is a forward stochastic differential equation, the SDE is 



Xl'^ = x + j^ dWt, t<s. 
3.4 Example and Simulation 

The structure of solution is interesting. Note that the collection {Tt, t E [0, T]} 
is neither increasing nor decreasing, and it does not constitute filtration. 

Example 3.3. 

~dYt = ZtdBt-ZtdWt, 
Yt = Wt- 

The solution is 

Yt = {BT-Bt) + Wt. 
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We usually apply binomial tree model to simulate Brownian motion. Wt is 
a forward binomial tree and Bt — i?t is a backward binomial tree. Then the 
coupled binomial tree is a tetrahedron. It is could be illustrated that all the 
paths {t, Wt, Bt — Bt) are within a tetrahedron. 

(Wt, Bt — Bt) is a coupled Brownian motion. 

Figure (P) illustrates Wt, Figure ^ illustrates Bt — Bt, and Figure ^ 
illustrates {Wt, BT — Bt). The tetrahedron is big and the paths arc concentrated 
by central limit theorem. 

4 Main Results: Convergence Results for Dis- 
crete BDSDEs 

4.1 Convergence of The Solution for Discrete BDSDEs 

We consider the discrete terminal condition is := ^" = <i>((WJ'), < j < n), 
which is CJ,"„-measurable random variable, for the discrete case. Firstly, for the 
scheme ^ of BDSDE, if wc construct the processes: 

then the convergence between (?/", z") to (yt, zt) can be derived in the same way 
as Donsker-Type theorem for BSDEs , by (P.Briand, B. Dclyon and J. Mcmin. 
(2001)@), 

Assumption {H.2) ^ is J^f-mcasurable and, for all n, ^" is t/^*„ -measurable 

s.t. 

£[^2] +supE[(rf ] < oo 

n 

Assumption {H.3) ^" converges to ^ in as n ^ oo. 

Theorem 4.1. If the assumptions (H.l), {H.2) and {H.3) hold. Let us consider 
the scaled random walks i?„, Wn, if B„ —i- B, Wn — i- W as n oo in the sense 
of that 

lim sup \Bt - =0 m P, 

"^°°0<t<T 

and 

lim sup \Wt-Wl'\^0 inP, 

n-*cxD o<t<T 

then we have (y", z") — > (y, z) in the following sense: 

lim I sup \yl' - ytl V / \zl' - z,|^ ds I m P. (12) 

n^oo I o<t<T J J 

Method for the proof. The key point is to use the following decomposition 

y" - r = (y" - r"-") + (r"-" - y°°-p) + {y°°^p - y), (13) 
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Figure 1: (t, Wt) 




Z" - Z = {Z" - Z'^'P) + - Z°°^P) + {Z°°^P - Z), (14) 

where the superscript p stands for the approximation of the solution to the 
BDSDE via the Picard method. More precisely, we set = 0, = 

0, F"'° = 0, Z"'" = and define (y°°^p+^ ^ z°°'P+^) as the solution of the BDSDE 

T T T 

y^oo,p+i =^+J fiY^'P, Zf'P)ds+ j .g(y,°°'P, Z^-P)dBs-J Z^^P+^dWs, 0<t<T. 
t t t 

(15) 

(•,^yoo,p+i^ ^oo,p+i) jg solution of a BDSDE with non-dependent but random 
coefficients) and similarly 

VV^^ = Vk+i^ + fiVk^^ 4^^)^ + aiVk+i, 4+i)'^^k+i - z'^-P+'^VsPk+i, k = n-l, 

n.n,p+l _ pn 
ifn S 

In order to define the discrete processes on [0, T] we set for < i < T, 
Y^'P = y]lfs] and Z^'P — z^^fg^ so that Y^'P is cadlag and Z"-'P caglad (cadlag 
means right continuous with left limits and caglad means left continuous with 
right limits). 

We shah prove in LemmalT^that the convergence of (F"'P, Z"'P) to (y", Z") 
is uniform in n for the classical norm used for BDSDEs which is stronger than 
the convergence in the sense of 1)12^ : this part is standard manipulations. 

We shall prove that for any p, the convergence of (F"-?', Z"-p) to (y°°^P, Z°°'P) 
holds in the sense of ([12^: this is the difficult part of the proof, and we shall 
need the results of section 14.1.11 

4.1.1 Convergence of Filtrations 

Let us consider a sequence of cadlag processes W"" = (Wj")o<t<T and W ~ 
iWt)o<t<T a Brownian motion, all defined on the same probability space (fi, P); 
T is finite. We denote by (C/") (resp. {Gt)) the right continuous filtration s.t. 
a{W") C CJ"(resp.CT(M^) C Gt)- Let us consider finally a sequence X" of Gj-' 
measurable integrable random variables, and X an t/x-measurable integrable 
random variable, together with the cadlag martingales 

M," =E[X" I g,"], Mt=E[X\Gt] 

We denote by [M", M"] (resp. [M, A/]) the quadratic variation of A/" (resp. 
M) and by [M",W"] (resp. [Af,l^]) the cross variation of A/" and W" (resp. 
M and W). 

Theorem 4.2. Let us consider the following assumptions 

(Al) for each n, is a square integrable G" -martingale with independent 
increments; 
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(A2) W" —^Win probability for the topology of uniform convergence of 
cadlag processes indexed by t € [0,T]; 
(A3) a. E|X2| +supE|(X")2| < oo 

n 

b. E[\X" -X\] 0; 
Then, if conditions (Al) to (A3) are satisfied, we get 

, M" , [Af" , M"] , [M" ^ {W,M, [M, M] , [M, W]) m P 

for the topology of uniform convergence on [0, T]. Moreover, for each t e [0, T], 
for each < S < 1, 

{Wr. , [A/", [Af", W"]]/') ^ (m, Af„ [A/, M]]/', [A/, W]'/^) m L'+\n,T, P). 



Corollary 4.1. Let W and VF", n G N*, be the standard Brownian motion 
and the random walks of Theorem 14.11 Let us consider, on the same space, X 
and X" satisfying the assumption (A3) of Theorem 14.21 

Then there exists a sequence (Z")o<t<T of CJ."-progressively measurable pro- 
cesses, and an fj. -progressively measurable process (Zt)o<t<T such that: for all 

te[o,T], 



t 



A/t" = E[X"] + j Zl'dW^, Mt = E[X] + j Z.dWs 





and 



T 



j (Zf - Ztfds 



Moreover, if < (5 < 1, Z" converges to Z in the space L'^+^{n x [0, T], : 
;B([0,r]),P®A) where A denotes the Lebesgue measure on ([0, r],S([0, T])). 

4.1.2 Proof of Theorem ICT 

Equations (fT3)) . ([T4)) with the following lemma proved in appendix A. 



Lemma 4.3. Here we need to assume that 



hm sup \Bt - -Bj"| =0 in P. 



With the notations following (fTS]) . ([16 



supE[ sup - Fj"-''!' + / |Z; - ds] ^ O, as p ^ oo 

n 0<t<T J 
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imply that it remains to prove the convergence to zero of the process y"'P — 
y°°'Pand — Z^'P. This will be done by induction on p. For sake of clarity, 
we drop the superscript p,set the time in subscript and write everything in 
continuous time, so that ([T5|) . (fTC|) become 

T T T 

y/ = e + J /(n, Zs)ds + J giY,, Z,)dBs - j ZidWs, 0<t<T. 



1 1 1 

Vr =C + J /(n- , ^Drf^s + / 5(n" , Z^)dB^ - J ZTdW^ , < t < T. 
t t t 

where A" = [s/^li^ and Y- denotes the caglad process associated to Y . The 
assumption is that {Y^ , ZJ^}o<t<T converges to {Yt, Zt}o<t<T in sense of 
and we have to prove that {F/", ^^"}o<t<T converges to {Y"/, Zj.}o<t<T in the 
same sense. 



According to the Peng and Pardoux's paper [8|], we define the filtration 
{Gt)o<t<T by 

Gt = V 
and the CJt-square integrable martingale 



T T 



Mt = E^* [i + j f{Ys, Zs)ds + j g{Ys, Zs)dB,], 0<t<T. 



Then there exists C/t-progressively measurable process {Z^} such that 

T 

E J \Z'fdt < oo 



t 

Mt = Mo + / Z'JWs, < t < T. 







On the other hand, the process, defined by 
t t 

+ J fiY:i , Zl^}dA^: + J g{Yr_ , Zl^)dBl^ , < t < T, (17) 



satisfies 



z 

Mj" = MS + J Z'.'^dW^. (18) 
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Hence M" is an jr."-niartingalc and, since = 

T T 

Mr = I gn, = vr + j f{Y^_,z^)dA: + J g{Y:i,zi^)dB^. 



(19) 

If we want to apply Corollary, we have to prove the convergence of M^. 
But since and Z" arc piccewisc constant, we have 



-Yt- JfiY,,Zs)ds - Jg{Y,, Z,)dB, 





< \Y^ ~Yt\ + J Z^) - /(n, Z,,)\ds + 



JgiYp, Z^)dB^ - ]g{Ys, Zs)dB, 





< {1 + KT) sup \Yt" -Yt\+K J\Z^ - Zs\ds^ 

0<t<T 



Jg{Y:\Zl^)d{B^-Bs) 



J{g{Yr,Z^)-g{Y,,Zs))dBs 





< (l + KT) sup \Y^!^ -Yt\+Kj\Zl' - Zs\ds + KVT sup 



0<t<T 



0<t<T 



J{K |y«| + a \Z^\ + |g(0, 0)|)d(i3," - S,) 



which tends to zero in probability and then in by L^-bounded. This and 
equations p8|) . (fTg]) . imply together with Corollary that M" converges to 



Mt=VP*[i+j f{Ys,Zs)ds+ j g{Y,,Z,)dBs\ - Y^+J f{Y,,Z,)ds+ J g{Y,,Z,)dB, 



in the sense that 

T 

sup |Aff-Mt|+ / |Zf -Z.l^ds^O in P. 

0<t<T J 


Since we want to prove that 

T 

sup IF/" -Y;\+ f\ZT-Z,fds^O in P, 

0<t<T J 



it remain only to demonstrate 



sup 

0<t<T 



r I r E 

1 f{Y:,Z:)dA:- J f{Y,,Z,)ds + J g{Y:_,Z:)dB:^ J g{Y,,Z,)dB, 



in P. 
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This is true since we have just proved the convergence of / — f{Ys,Zs 



to zero in probabihty and since the jumps oi t J f{Yp, Z]!;)dA^ tends to zero 



according to 

t t 



sup 

0<t<T 



in P, sup 

0<t<T 



J {z:rdA': - J zids 



in 



while we also have proved the convergence of 
to zero in probability. 



Sg{Yp,Z^)dB^-Sg{Y,,Z,,)dB, 



4.2 Convergence of Modified Solution 

Theorem 4.4. If the assumptions (HA), {H.2) and {H.3) hold. We also as- 
sume that 

lim sup \Bt -B]^\=0 m P. 

Then the discrete solutions {(y", z")}^i under the scheme Q converge to the 
solution {y, z) of ([7]) in the following senses: 



lim <^ sup \y1-yt\'^+ / 

i^oo o<t<T J 



- ds ^ = mP. 



(20) 



This can be derived from the convergence (|12|1 . 

For the convergence of this scheme, we must consider the following estimates 
under the following: 

Assumption C" G L^iC), IE[E 1/(0, 0)|'] < oo, |.9(0,0)|'] < oo. 

3=0 ]=a 
For this reason, we need the following Gronwall type lemma, which is proved 
in 3. 

Lemma 4.5. Let us consider a,b,a positive constant, M < 1 and a sequence 
(fj)j=i^...^„ of positive numbers such that, for every k 



then 



Vj + a < a + bS Vi , 



supvj + a < aes{b), 



(21) 
(22) 
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where es(b) is the convergent sequence: 

°° hP 

ssib) = l + Y,-{l + S)---{l + {p-l)S) (23) 
p=i P- 

which is decreasing in 6 and tends to as J ^ oo. 

Lemma 4.6. We assume that 5 is small enough such that {1 + 2K + 7K'^)d < 1. 
Then 

n 

supE ly^l' + 6E[Y^ |z;f ] < Q",/,ge(i+2^+^^')^ (24) 

^ 3=0 

where = |/(0, 0)|^ + 3 |5(0, 0)|^ + (1 + K6 + 3K^6 + 3a^VS)E ■ 

Proof. By explicit scheme, we have 

We then have 

E\y]\^-E\y]^,\^ = -E\z'l\'s + E[giy]+„z^^,)]H-E\fiE[y]+, \ Q^.lz])]' 6^ 
+2E[y^-f{E[y'^^,\g^^],z])]6. 

(25) 

Taking sum for J = z, n — 1 yields 

71 — 1 71—1 

E < E\a^-Y,E\z]\'s + 25EY,{\y]\ (1/(0, 0)| + K \E[y]+, \ G^^] \ + K \z^\)} 

n-1 

+SEY,MOM + K\y]+^\+a\z]+,\)\ 



Since the second last term is dominated by 

n-l 



n— 1 ^ 

,5E^{|y7|' (1 + X + + 1/(0,0)1^ + A' \E[y-^, \ g-]\' + - \z]\'} 



3=1 



n-^ 

< <5E ^{|y;f {1 + K + 4K^) + 1/(0, 0)1 V - l^f} + KSE \e\ 

3=i 

and the last term is dominated by 

TI-l 

3<5E5^(|5(0,0)|Vi^^ + 



3=1 
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< 3<5E^(|g(0, 0)1' + \y^\' + jz^"!') + S^K^S + a^V6)E \C\ 



3=1 

we thus have 

n-l 



E|yj"|'+5(^-3a2)^E|z7|' < \f {0,0)f + 3\g{0,0)f + {I + KS + 3K^S 



n-l 

2 



y 



Then by Lemma I4.5i we obtain ([24]) . 

Proof of Theorem 14.41 The convergence of {y"',z") to {y,z) is aheady 
proved above. To prove that of (y", z"), it is sufhcient to prove 



z'^-z^l^ds ^ =0. 



hm J sup iyr-2/rr+ / 
y. 

From ([6]) and we have 

E I - 2/;.f = E I y^Vi - I ' - -^E \z] - | ' + <5E | giy]+, , z^Vi ) - givj+i, z^+i) 

-E 1/(2/;, z;) - /(E[y^i I z;Or -5' 

+2E[(y;' - ■ (/(y^, z^) - f{E[y-^, \ g^,lz]))]S. 

We then take sum over i from j' to n — 1. With — £, = 0, we have 

n — 1 n — 1 

E|yr-yri' < -5E^|z;-z7|' + <5^E|5(f^i,z7+i)-5(y;+i,z;+i)|' 
+2<5 ^ E[(f7 - y-) ■ {f{y-, z^) f{E[y^^, \ G"^,],^))] 

n — 1 71—1 n — 1 

< -SE 1^7 - 41' + 2A'25 ^ E IfJ+i - y^+i f + 2aH ^ E [z^+i - 

n—1 71 — 1 

+2if25 ^ E I - 2/7 1 V (5/2 ^ E I 



I z" - z" 



n-l 
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Since y^~E[y-_,, \ G^^] = f{E[y^^, \ g^^lz^)5 + VSE[g{y]^„z^^,) \ 
^"j]£j+i, the last term is dominated by 



n-l 

|2 



j = i j=i 

n-l 

+ ^E|E[g(f;Vi:^"+i)l^,"]r^- 



3=1 



But with (|24p. the second term is bounded by CS^, and the last term is 
bounded by 



n—l n— 1 



4K^S J2^\y]\ + 4a''5 ^ E \z]\ + 4{K^S + a^V5)E |f f + 2 |g(0, 0)|' . 



We thus have 

n—l n—l 

'2 



^ n—l n—l 

E m - y"f + h- 2a^)SE J2 M - 41' < i^K^S + K^S^ + 5^ + 2KS^) ^ E - y^j 



n—l n—l 



+4i^2 J ^ E I I + 4a2,5 ^ E I I 

+4(^^25 + a"V^)E iri' + 2 |g(0, 0)|' + C-J^. 
According to Lemma 14.61 and here providing that g(0, 0) = 0, we further have 

n—l n—l 

E lir - + (2 - 2a')'5E^ |z^ - z,"!' < (1 + 2fc + Sfc^)^^ E j^;^ - y^f + C'VS. 

j=i j = i 

By Gronwall Lemma 4.2.2, we get 

supE - < C'V^e(i+2'=+5'=')^. 

i<n 

Then these two inequalities implies (|20|) . 

Appendix A. Proof of Lemma 14.31 

For the proof of this lemma we come back to the discrete notations and we show 
that 

Lemma A.l There exist a > 1 and no £ N such that for all n > hq, for all 

^3 ^ 
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where, for p G N, 



sup a 

0<k<n 



kS 



Vk -Vk 



Tl-1 

k=0 



kS 



^k ^k 



Proof. For notational convenience, let us write y, z in place of y 



Ji,P+l _ ytl-P 



^n,p+i _ ^n,p g^^^ y place of y^'P — y^'P~^, z^'P — z"'P~^. Let us pick ip > 1 
to be chosen later. With these notations in hands, we have, for k ~ Q, ?? — 1, 
since = 0, 



We write - yf^^ = 2yi{yi - yi+i) - {y, - yi+if , to use the equation (HI]), 



since 



y^-y^+l = 5{/(?;f''',zf'P) - /(j/f'P \zf'^ ^)} + ^{9iy?+i, - g{y"+^i \ ^)}ej+i 



According to (H.l), we have, for each u > 0, 

< 2iKy,y)yf + ,yiuj+vf), 

'^Viigiyifi, z-fi) - 5(yr/i~\ < 2 ly^l (iC + a \vi+i\) 

< [{K' + a')/,^]y^ + u{u^^,+vf+,) 

and moreover, (A.l) implies easily that 

Szf < 3(y, - y,+if + &KH\u^ + vf) - 6S{K^u'^^, + a^v^,). 
As a byproduct of these inequalities, we deduce that, for k ~ 0, n — I, 



(A.l) 



i—k 

n-1 

< 2K'{6/:.) ^ + E ^'("? + + [^(^' + "')/'^] E ^'y's^+i 

n— 1 n— 1 



n-l 



i—k 



n-l 

i—k 



n-l 



n-l 



n-l 



i—k 



i—k 



i—k 
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and, setting p = {v + 2K'^S)ipS, we get 

n-l 



71-1 



n-l 



i—k 



■ ^{5/3) ^^zf < (1 - ^ + 2K'6^/^) J2 f'yf + P E + ^'•') 

i—k i—k 
n — 1 n — 1 

+[V5^{K^ + c?)lv\ E (^'y- + i^<^V^E ^'("^+1 + 

i— A; i—k 
n— 1 n— 1 



i—k 



Thus, if 1 - + 2K'^6iplv < 0, we have, for fc = 0, n - 1, 

n-l 

i—k 

n—1 n—1 

i—k i—k 

n~l n—1 

i—k i—k 



(A.2) 



In particular, taking the expectation of the previous inequahty for fc = 0, we 



get 



n-l 



n-l 



1=0 4=0 

Now, coming back to {A.2), we have, since j/„ = 0, 

n-l 

sup < pE^'("?+^^') + 2[^</'(^^' + a')/H_fup 

n-l 

Y,^'iu^+l+V^+l)^^+l 



(A.3) 



0<fc<n 



i=0 



0<*;<n-l 



n-l 



J2 '^'yh'^+1 



-2vLpV5 sup 

0<fe<n-l 



1=0 



'iip^fS sup 

0<fc<n-l 



n-l 



and using Burkholdcr-Davis-Gundy inequahty, we obtain, 

n-l 

0<A;<n 



n-l 



n—1 Ti— 1 

4=0 i=0 

n-l 

< pE[E V'X^i' + + Ci^[(/^2 + a2)/HE[ sup 



0<A:<n 
n-l 



+C2;^<5E[E^'(^m+^^4+i)] + C3V'5IE[Ev''^'] + ^E[ sup ip^yi] 



4=0 



4 = 



0</£<T 
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Finally, from {A.3), we get the inequality, 



n—1 n—1 



E[ sup + J ^ < AE[ sup ^^'^.^ + <5 ^^z;,^]. 



0<A:<n 



i=0 



0<A:<n 



1=0 



y(i/+2X^(5)(l+3C|y-3Ci(-ft:''+a^)/i/)+C2i^ 



(A.4) 

and providing that 1 — <^ + 



where A 
2K'^5ip/v < 0. 

Firstly, we choose ly such that '^(i+'^'^3~3C^i(-^ +a )/v)+C2i' _ ^12^ We con- 

'^—Ci(K-'+a-')/i' ' 

sider only n greater than ni (i.e. K5 < 1 and 2K^5/v < 1). Let us pick 
of the form 7* with 7 > 1. We want that 1 — 7'' + 2K'^5^^ /v < meaning 
that 7 > exp-5-^log{l - 2K'^6/v). Since exp{-6-Hog{l - 2K'^6/v)} tends to 



exp{2K'^5/v} as 



J ((5 0), we choose 7 = exp{l + 2K^5/v}. Hence, for ' 



greater than 712 the condition is satisfied and (4.18) holds for (/3 = 7"^. It remains 
to observe that, v and 7 being fixed as explained above, A converges, as n — > 00, 

which is equal to 1/2. It follows that for n large 



to 



y(l+?,C^-?,Ci{K^+a^)/v)+C2y 
l-Ci{K^+a^)/v 



enough, say n> no, X < 2/3 and 



n-l 



n-1 



E[ sup 7''y? + S Y'zf] < 2/3E[ sup -f^'ul + 5 ^ 7^^^ 



0<fc<n 



i=0 



0<fc<n 



i=0 



which concludes the proof of this technical lemma. 

To complete the proof of Lemma 14. 3[ it remains to check that 



supE 



I n.l|2 

sup I 

0<fc<n-l 



i=0 



is finite. But it is plain to check (using the same computations as above) that 
for n large enough. 



E 



I n.l|2 

sup \y^' I 

0<fc<n-l 



i=0 



<Ce'- (E[e2] + |/(0,0)p + |5(0,0)|2) 



where C is a universal constant. 
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